The transition to chaos of coupled oscillators: An operator fidelity susceptibility study 
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The operator fidelity is a measure of the information-theoretic distinguishability between perturbed and unper- 
turbed evolutions. The response of this measure to the perturbation may be formulated in terms of the operator 
fidelity susceptibility (OFS), a quantity which has been used to investigate the parameter spaces of quantum 
systems in order to discriminate their regular and chaotic regimes. In this work we numerically study the OFS 
for a pair of non-linearly coupled two-dimensional harmonic oscillators, a model which is equivalent to that of 
a hydrogen atom in a uniform external magnetic field. We show how the two terms of the OFS, being linked 
to the main properties that differentiate regular from chaotic behavior, allow for the detection of this model's 
transition between the two regimes. In addition, we find that the parameter interval where perturbation theory 
applies is delimited from above by a local minimum of one of the analyzed terms. 
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I. INTRODUCTION 

Given a classical system in which chaotic dynamics arise, a 
fundamental problem is to understand the behavior of the cor- 
responding quantum analogue of this system. A starting point 
for this comparison is the observation that while classical dy- 
namical chaos can typically be described in terms of the diver- 
gence of initially neighboring trajectories in phase space, the 
unitary evolution of closed quantum systems does not allow 
such a characterization. So, how does classical chaos mani- 
fest in a quantum system? Is quantum chaos reflected in the 
distribution of the energy levels of the corresponding quantum 
system, for example, or in terms of the temporal evolution of 
suitable expectation values? For the last several decades, in- 
tense effort has been devoted to the study of quantum chaos 
111]. A variety of approaches have been taken, including ran- 
dom matrix theory [2], quantum motion reversibility [3], sta- 
bility [4], fidelity [5], entanglement [6], and recently measures 
of phase-space growth rates nfl. All of these techniques suc- 
cessfully address specific aspects of quantum chaotic behav- 
ior. The fidelity, a quantity commonly encountered in the con- 
text of quantum information theory, has been used extensively 
in the form of the Loschmidt echo to study quantum chaos by 
Prosen, et. al [5] and others. See Ref. [8] for an extensive re- 
view of the fidelity approach to quantum chaos. The response 
of this quantity to infinitesimal perturbations, what we call the 
operator fidelity susceptibility (OFS), was formulated from a 
differential-geometric perspective in ifloL fllfl - This quantity 
has been studied in the context of quantum chaos in [9], where 
the state- and level-dependent contributions to the OFS were 
distinguished. The operator fidelity and OFS are the gener- 
alizations of the ground state fidelity and fidelity susceptibil- 
ity, resepectively, that have been fruitfully used [12-14] over 
the last years to characterize another important class of phe- 
nomena encountered in many-body quantum physics: quan- 
tum phase transitions. The state fidelity approach is based 
on the idea that one can detect quantum critical points by 



employing a geometric measure of distinguishability between 
ground states corresponding to neighboring parameters. The 
OFS is a natural extension of this approach from the space 
of states to operators. In an analogous fashion to the ground 
state case, the OFS gives a measure of the rate of distin- 
guishability between "neighboring" members of a family of 
operators. Here we consider the unitary evolution operators 



Ux 



-iH(X)t 



generated by a family of time-independent 
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Hamiltonians H(X) = H + XV. With this choice, the OFS 
measures the rate of separation between the unitary evolutions 
induced by U\ and Ux+sx, where 5 XV is an infinitesimally 
small perturbation to the Hamiltonian. 

The reason why the OFS xW allows one to distinguish 
the transition from a regular to a chaotic regime of a quan- 
tum system is two-fold. First of all, the theoretical study of 
x(A) carried out in I9D and based on random matrix theory 
arguments [2] has shown on general grounds that for systems 
characterized by random perturbations V drawn from proper 
ensembles, the quantum chaotic evolutions can be character- 
ized as those which have the highest resilience to these pertur- 
bations. This can be seen by expressing the OFS in terms of an 
autocorrelation function of the perturbation |8]. In this form, 
when the system is chaotic the correlation function may decay 
more rapidly than in the regular regime, leading to a slower 
decay of the fidelity. Furthermore, as described in [91 lllll . 
and as will be reviewed in the next section, the OFS can be 
split into two terms, x(A) = X (ty + X (A), one which 
depends on the variation of the energy levels and the other on 
the variation of the eigenvectors. In particular, the first term 
depends on the variation of the eigenvector and is directly re- 
lated to the spacings between energy levels. In this respect, 
y^fA) is analogous to the ground-state fidelity susceptibility 
ltl2l[l4ll . a quantity which depends on the energy gap between 
the ground and first excited states. This property links the 
OFS to one of the main paradigms for the definition of quan- 
tum chaos: the statistics of the spacings s n = E n+ i — E n 
between neighboring energy levels. Indeed, in a seminal pa- 
per Bohigas, Giannoni, and Schmidt |15] proposed that cer- 
tain universal level spacing statistics encountered in random 
matrix theory ||2|] also arise in the spectra of regular or chaotic 
quantum systems. In particular, regular systems are expected 



to have Poissonian level spacing statistics, P(s) = e~ s . In 
this case the most probable value is s = 0, since frequent 
crossings between members of uncorrected subsets of energy 
levels may occur. This property is a direct consequence of the 
integrability of the system, corresponding to the existence of 
a sufficiently large number of conserved quantities. On the 
other hand, with the exact form of the distribution depending 
on the few existing symmetries of the model, chaotic systems 
may display Wigner-type statistics, P(s) sa s^exp— s 2 . In 
this case level crossings are suppressed as a result of correla- 
tions between the energy levels due to a lack of symmetry in 
the system. The sensitivity of the eigenvector part of the OFS 
to these level statistics has already been successfully tested 
in 13], where it described the transition to chaos in the Dicke 
model. The same basic ideas have subsequently been used to 
define a simplified version of the OFS that has been applied 
to spot the avoided crossings in some Bose-Hubbard systems 

In this paper we employ the operator fidelity susceptibil- 
ity to investigate a system of two non-linearly coupled two- 
dimensional harmonic oscillators. The importance of this 
model arises from its equivalence to a prototypical system 
for which the regimes of classical and quantum chaos have 
been well-documented theoretically, and for which experi- 
mental data exist and agree exceptionally well with theoretical 
predictions: the hydrogen (or hydrogen-like) atom in a time- 
independent and uniform external magnetic field lfl7l - [l9ll . In 
order to spot the transition from regularity to chaos, we com- 
pute for the oscillator problem both parts of x(A). Our analy- 
sis, besides providing a further test of the applicability of the 
OFS, will allow us to compare their different behavior in the 
transition and to give an interpretation of the full OFS. 

In Sec. |II] we briefly review the operator fidelity and op- 
erator fidelity susceptibility. In Sec. [TTT] we explain how the 
model is derived, describe the algebraic representation allow- 
ing for its efficient diagonalization, and detail the numerical 
procedure adopted. In Sec. [IV] we outline the level spacing 
statistics for the model and give the results of the operator fi- 
delity susceptibility analysis. We conclude with Sec. [V] 



II. OPERATOR FIDELITY 

Due to the preservation of inner products under unitary evo- 
lution, the classically useful notion of the divergence of tra- 
jectories in phase space resulting from sensitivity to initial 
conditions does not apply for quantum systems. One may 
instead compare the unitary evolution operators correspond- 
ing to nearby points in parameter space. This approach al- 
lows the problem of discriminating between the regular and 
chaotic regimes of a given quantum system to be cast into the 
information-theoretic idea of statistical distinguishability be- 
tween states representing evolution operators. 

We start our discussion by recalling that the overlap be- 
tween quantum states belonging to a /^-dimensional Hilbert 
space H generalizes to the space C(H) of linear operators 
acting on H in a simple way. Supposing, for example, that 
p = J^i Pi\i)(A i s m e density matrix of a given quantum sys- 



tem, one can identify any operator X G H with the state 
\X) e T-L (£> % defined in terms of the purification of p, 



| A) = A®1\ 



(1) 



The p-fidelity between two operators X, Y is then the (bi- 
partite) state fidelity: 



F p {X,Y) = \{X\Y) p \ = \Tr{pX^Y)\. 



(2) 



where (X\Y) p defines an hermitean scalar product over 
£("H). In this work, as in J3, [HI], our starting point is to 
compute the p-fidelity between members of the one -parameter 
family of unitary time-evolution operators, U\(t) = e~ ltHx , 
corresponding to slightly different Hamiltonian parameters 
A, A + SX. The second-order expansion \U\+sx)p = \U\) p + 
SX\d\U\) p + 8\ 2 \d\U\) p /2 of the state representing the per- 
turbed evolution, together with the relations obtainable by the 
identity d\(U\\U\) p = 0, allow the second-order expansion 
of the fidelity to be written 

SX 2 
F P {Ux, Ux+sx) = 1 - — x P (A), (3) 

in terms of the operator fidelity susceptibility (OFS) llXltl - 

Xp (X) := {d x U\dxU) p - \{d x Ux\Ux) p \ 2 . (4) 

In this work we always choose states p which commute with 
the Hamiltonian, Hx = Hq+XV. The spectral decomposition 
of p in a basis of energy eigenstates of Hx and the application 
of time-dependent perturbation theory 111111 allows the OFS to 
be decomposed into two parts, x p (A) = \p (A) + Xp W> 
where [20] 

D 

X «(A) = 27rf Y, Pn\(n\dxH\m)\ 2 G t (E n -E m )(5) 



D 



X P 2) W = t 2 (J2pn\(n\dxH\n)\ 2 
- \J2Pn(n\dxH\n)\ 2 ), 



(6) 



D is the dimension of the Hilbert space, and Gt(x) = 
sln 7T ^ x i ■ Following the discussion in Ref. lHHl . Xp is 

(2) 

associated with the change of the eigenvectors and xh m e 
eigenvalues [20]. Note that the OFS may alternatively be 
expressed in terms of a dynamical two-point autocorrelation 
function of the perturbation 1 8] . 

(1 2) 

From the explicit form of xp' (A), one may identify how 
the OFS reflects the chaotic behavior of the system. In par- 
ticular, we will be interested in the long-time behavior of the 
OFS. Observe that for xp (A), at large times t the function 
Gt{x) acts as a filter selecting the contributions due to neigh- 



boring levels in the sum ©. In fact, since 



fitix) 



5(x), the largest contributions to Xp come from small gaps 



(i) 



s n = E n+ i — E n < l/t. Therefore, for large t Xp W be- 
comes highly sensitive to the specific level spacing distribu- 
tion that characterizes the model for a given A. As for the 
second term of the OFS, it displays a quadratic dependence 



unstable. For the classical model, it can be seen from scal- 
ing arguments that the degree of regularity or chaos is deter- 
mined chiefly by the single parameter ? — E^ 2 ^ given by 
the energy E and magnetic field strength 7, where the limit- 
ing cases are e — > — 00 (e — > 00) for the Coulomb (Landau) 
regimes 11711 . The degree of chaos is dependent on the rel- 
ative strengths of the Coulomb attraction to the nucleus and 
the diamagnetic interaction of the electron with the magnetic 
field. For e less and less negative, the perturbation due to the 
magnetic field begins to dominate, and the phase space be- 
comes increasingly chaotic. Note that the analysis here con- 
siders only the bound spectrum. 

Returning to the quantum mechanical hydrogen atom, with 
tor G t (E n - E m ) in X p ' have a narrow-enough peak that it a suitable change to " se miparabolic" coordinates n = j7T~z 



on time and may be written as xp (ty — t 2 A 2 Vd, where 
A 2 Vd = (Vd 2 ) p — (Vd) 2 , is the variance of the diagonal part 
Vd of the perturbation V in the energy eigenbasis. 

We remark here on the domain of validity of the OFS ap- 
proximation for the fidelity. In order for the expansion (O to 
be valid, it is necessary that T p ~ 1. Since the rate of growth 
of the OFS is at most quadratic in time, it is thus necessary that 
(tSX) 2 <C 1. Provided this holds, it is also desired that the fac- 

,(!) 



samples primarily the nearest-neighbor level spacings. If this 
is satisfied, xp is expected to be sensitive to the level spacing 
statistics. Given the set of energy levels in the support of p and 



defining some typical level spacing A 



iyp 



(E 



n+l 



E, 



n/typi 



the time scale in which Gt (x) samples primarily the nearest- 
neighbor level spacings is given by l/t w A typ . 

In our analysis we consider two kinds of initial states, re- 
stricting ourselves to a single parity sector of the model in 
order to allow comparison with the universal level spacing 
statistics. On one hand we will use a maximally mixed state 
po = i/D, defined as a truncated version of the exact state p 
whose range will be the £>-dimensional space of numerically 
well-converged states in the even sector. On the other hand 
we consider p to be a (truncated) Gibbs thermal state over the 
even sector, p\ = e~^ Hx /Z\, where Z\ = Tx(e~^ Hx ). This 
will allow us to establish the extent to which the introduction 
of temperature in the system modifies the behavior of the OFS . 



III. THE MODEL 



and v = y/r — z, the Schrodinger equation for the hydrogen 
atom of energy E may be written 



[^+^-^ 2 +^) + ^M 2 ^ 2 (M 2 + ^ 2 )-2]|*) =0, (8) 



where D,, 



i- ~Hw? + jw)' and similarl y for D »- 

If one now performs another coordinate transformation to 
"dilated" semiparabolic coordinates ilal . u = (— 2_E) 1 / 4 /i 
and v = (—2E) 1 / i v, the Schrodinger equation takes the form 



[D u +D v - 



7 



8(-2E)' 



2 21 2 1 2\ 

u v (u +v )- 



V^2E 



]|*) = 0. 



Here it can be seen that the hydrogen atom is equivalent to 
a pair of coupled two-dimensional harmonic oscillators, each 
having identical angular momenta which are in this case zero. 
If we now define A = -f 2 /(-2E) 2 and e = (-2E)- 1 / 2 , the 
final form for the Schrodinger equation is 



The hydrogen atom in a uniform external magnetic field 
is one of the simplest time-independent systems exhibiting 
quantum chaos, and has been studied in detail analytically, 
numerically, and experimentally [ 17-19]. The Hamiltonian of 
a non-relativistic hydrogen atom in a uniform cylindrically- 
symmetric external magnetic field has the form 



rA, 



H= Y -- + -L Z + Y ( X +y), 



7V2 



(7) 



where 7 is a dimensionless parameterization of the magnetic 
field strength and L z is the component of angular momentum 
along the magnetic field axis. This component of angular mo- 
mentum is conserved, and we restrict ourselves to the L z = 
sector. For weak magnetic fields a non-zero angular momen- 
tum quantum number simply results in a uniform shift of the 
energy levels within a given sector, the Zeeman effect. 

The application of a magnetic field to the classical hydro- 
gen atom gives rise to the possibility of richly complex elec- 
tron orbits. Two principal orbital modes of the electron can 
be seen: a "rotational" mode in the z = plane and a "vibra- 
tional" mode along the z-axis [17]. An effect of a sufficiently 
strong magnetic field is to make trajectories near these modes 



D r 



<^-^ + ±u 2 v 2 (u 2 +v 2 ) |*>=e|*>. (9) 
4 16 J 

The problem we aim to solve in this work is to find the 
eigenvalues e corresponding to the choice of coupling param- 
eter A. This is called the oscillator problem. Notice that our 
choice of A fixes the ratio j/E, so the oscillator problem es- 
sentially corresponds to finding the hydrogen atom energies E 
intersecting a curve of constant j/E in the j-E plane. Alter- 
natively, one may solve the hydrogen problem for the energies 
E for a given magnetic field 7. This is a generalized eigen- 
value problem of the form A\ty) = aB|*) (Eq. [HJ. In Ref. 
01811 both approaches have been taken, with similar results for 
the level spacing statistics. 



A. Algebraic representation and numerical procedure 

The problem of finding an efficient solution to the eigen- 
value problem for the hydrogen atom and to the equivalent 
nonlinearly coupled harmonic oscillators has been discussed 
in several papers. In particular, Delande and Gay were able 
to solve the problem elegantly by giving a dynamical group 



approach representation of the model IU9I1 . In the following 
we briefly review their results. One begins by defining the 
operators 



1 r , < & 1 d \ 
4 ou z u ou 

2 L <9u J 



(10) 



that generate the so(2,l) algebra and fulfill the commutation 
relations IU9I1 

(11) 
(12) 

(13) 

In terms of these operators, the Schrodinger equation (O can 
be written as 

A. 



[&xi &y\ 


= —iS; 


[Sy,S z ] 


= lb x 


[S Z ,S X ] 


— iby. 



[S z +T z + -V-e]\y)=Q, 



where 



V = (S x + S z f(T x + T z ) + (S x + S Z )(T X + T z f 

and T X:V , Z has the same form as S x , y ,z> except with u replaced 
by v. Note that [T a , Sp] = V a, e {x, y, z}. With this 
representation, the natural basis to use is the tensor product 
of the eigenbases of S z for each oscillator. That is, the set of 

states {\n} <g> \m)} n ,meN where 



Sz\n) = (n+-)|n) (14) 

S+\n) = (n + l)|n + l) 

S~\n) = n\n — 1) 

and S f± = S x ± iS^. In this representation, 5 Q — > 5 Q I 
and T Q — > I<8> 5 Q . The Hamiltonian for the oscillator problem 
may then be written 



where 



H x = S z ® 1 + 1 ® 5, + AV, 



V = - [E 2 (g) E + E <g> E 2 ] 



(15) 



(16) 



and E = 5 K + S z = S z + \{S+ + S~). Notice that the 
Hamiltonian is symmetric under interchange of the subsys- 
tems. Namely, defining the parity (or swap) operator P acting 
in the oscillator occupation basis as P\n) <E> \m) — \m) ® |n), 
we have [H, P] = 0. 

With this in mind, we now define an orthonormal basis 

U ,_[ MW)®\m) + \m)®\n)) {n^m) 



\n) ® \n) 



(n = m) 



\d n , m ) = — =(\n) ® \m) - \m) ® |n)) (n^m), 
V2 



where {|e nm )} ({|<in,m,)}) span the even (odd) parity 
subspaces. Since the parity operator P and Hamiltonian H 
commute, the Hamiltonian can be block diagonalized over 
the even and odd parity subspaces. From now on we restrict 
ourselves to states p which have support only over the even 
subspace, in order to better compare the fidelity analysis with 
the level spacing statistics. Indeed, in general for studies of 
quantum chaos it is useful to identify the symmetries of the 
problem and focus on individual symmetry subsectors. In the 
chaotic regime, for example, if the energy levels from several 
symmetry subspaces are included the level spacing statistics 
may take a non-universal form arising from the superposition 
of respectively universally-distributed subsets of levels. 

While in general the system to be analyzed has to be repre- 
sented in an infinite-dimensional Hilbert space H, in order to 
handle the problem numerically one must truncate the original 
Hamiltonian into one which acts on a finite-dimensional space 
Hk and appropriately represents the low-energy physics of 
the exact model. The above described algebraic representa- 
tion of the pair of two-dimensional harmonic oscillators easily 
allows for such a truncation. Here, the Hilbert space trunca- 
tion is implemented as an upper bound on the total allowed 
oscillator occupation, requiring n + rn < K, where K is our 
truncation parameter. The following is the truncated Hilbert 
space dimension for the even sector in terms of K: 



Dim(HT n ) 



'(f+ir 
m( z 2 



A'+l 



K even 
1), if odd. 



However, the eigenenergies of the truncated Hamiltonian 
poorly approximate the exact values near the top of the trun- 
cated spectrum. We address the inaccuracies introduced by 
the truncation by extending the truncation and measuring the 
corresponding changes in the energy spectrum, so that for a 
given truncation size we may determine the set of levels which 
are sufficiently well-converged to the true values for our pur- 
poses. We find that as the coupling grows the necessary trun- 
cation dimension also must increase to provide the same de- 
gree of convergence. This effect is due to shielding from the 
influence of the magnetic field due to strong coulomb attrac- 
tion for lower energy levels. For weaker magnetic fields a 
larger number of eigenvalues will be well-converged, but to 
be consistent we include only those levels which are well- 
converged over the entire coupling interval of interest. 

Recall that the parametrization of the coupling between the 
oscillators, A, is related to the hydrogen atom problem through 
A = 7 2 /(— 2E) 2 , where 7 and E were the external magnetic 
field strength and energy of the hydrogen atom, respectively. 
In this work, since we are studying the oscillator problem we 
only vary the oscillator coupling A and do not independently 
adjust the external magnetic field parametrizing the hydrogen 
problem. As previously mentioned, the quantity 77 = Ae 2 was 
found to characterize the degree to which the phase space of 
the classical model is chaotic [18]. In particular, for 77 < 1 
the phase space is almost entirely regular while for 77 sa 60 
the phase space becomes entirely chaotic [18]. It is known for 
this model that the character of the level spacing statistics is 
reflected in the proportion of the classical phase space which 



is chaotic. As we will observe now, where 77 is small the level 
statistics is Poissonian, P(s) = e~ s , while for larger 77 the 
level statistics transforms into a Wigner-Dyson form, P(s) = 



IV. RESULTS 
A. Level spacing statistics 

With these well-converged levels in hand, we now compute 
the level spacing statistics and compare them with the results 
of Delande and Gay [18]. Gradually increasing the coupling 
from the perturbative regime (Fig. QJ) through the so-called 
n-mixing regime, where levels from neighboring sectors be- 
gin to cross (Fig. [Tp) until the quantum chaotic regime (Fig. 
[TJ; andQJl), the level spacing statistics smoothly evolve from 
Poisson-like to Wigner-Dyson-like. Note that whereas the 
level statistics presented in [18] are the superimposed statis- 
tics for several coupling strengths, we have evaluated the level 
spacings for individual couplings. Our statistics include 800 
well-converged eigenvalues, with a truncation of K — 120. 
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FIG. 1: (Color online) Level spacing distributions for: (a) A = 
2.5 x 1CP 4 , approximately the location of the local minimum of 
X (1) (A), (b) A = 1.75 x 10~ 3 , approximately the peak of x (1) ( A )> 
(c) A = 3.0 x 10~ 3 , in the tail of the peak, and (d) A = 1.0 x 10 -2 , 
far into the quantum chaotic regime. Here all 800 well-converged 
levels are included in the statistics. The Poisson statistics are given 
by P(s) = e~ 3 (dashed curve) and the Wigner-Dyson statistics are 
P(s) = fse~^ s (solid curve). 



In Fig. |2]the eigenenergies are plotted as a function of cou- 
pling for two representative intervals of couplings and ener- 
gies. The changing character of the energy spectrum as a func- 
tion of the coupling strength is evident. For small couplings, 
degeneracies are lifted and levels from neighboring sectors be- 
gin to cross as the n-mixing regime is entered. As the coupling 
strength increases further still, the regular n-mixing regime 
evolves into quantum chaos as level avoidance begins to dom- 
inate. 



B. The x (1) term 

Having reviewed the behavior of the level spacings, let us 
now examine the operator fidelity susceptibility. We start by 
analyzing x > fi rst discussing the cut-off approximation: 



X {1) (Dc) = 



2lTt 



D, 



Z(D C ,(3) 



^2exp(-fiE n )C Do (n), (17) 



where C^ra) = E™= w „ \(n\V\m)\ 2 G t (E n -E m ),/3 = 
1/T, and Z{D C ,(3) = ^=1 exp(-0E n ). Here {\n}} are 
the (-y + l) 2 -dimensional or (^j-^)(^-^ + l)-dimensional 
(for K even, odd respectively) energy eigenvectors of the even 
parity block of the Hamiltonian Hk- The number of states in- 
cluded in the sum is called the cut-off dimension, D c . Note 
that if the odd sector were included as well, since the pertur- 
bation V preserves parity the sum in (fTTT i would be simply the 
linear combination of x computed over each parity sector. 
The dependence of all factors on both the coupling A and trun- 
cation dimension K has been left implicit. In our calculations 
we include only the states corresponding to well-converged 
eigenvalues, namely levels which vary little under additional 
growth of the truncation dimension K. Taking K = 120, for 
example, we retain D c = 800 out of the 3721 total levels in 
the even parity sector for the strongest coupling considered, 
A = 10~ 2 . Again, though the number of well-converged lev- 
els is maximal for smaller couplings, for consistency in our 
calculations we include the worst-case set of levels for the en- 
tire region of interest. 

For finite temperature, contributions to the OFS from higher 
energies are exponentially suppressed due to the coefficient 
p n = exp(—(3E n ) in the outer sum of Eq. ( TTTb . This serves 
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FIG. 2: (Color online) Eigenenergies vs. coupling strengths for K = 
120 for (a) couplings < A < 2 x 10" 3 and energies 30 < E < 34, 
such that < 77 < 2.3, and (b) couplings 5 x 10~ 3 < A < 10 -2 , 
such that 10 < 77 < 23. 
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(Color online) Partial sum 



X W {D,D) as a 
4.5 and A = 



function of D, for temperature T 

9.2 x 10" 5 , 2.5 x 10" 4 , 2.0 x 10" 3 , and 1.0 x 10' 



to give more weight to the low-energy part of the spectrum for 
low temperatures, while for high-enough temperatures all lev- 
els are included. In Fig. |3]we show the partial sum x^'(D c ) 
as a function of D c for up to 800 well-converged energy lev- 
els, for four different couplings A. For example, we can see 
that a temperature of about T = 4.5 is sufficient to provide 
a well-converged result for x*- 1 -', with the result being better- 
converged for larger coupling. This is due to an expansion of 
the energy spectrum towards larger energies as the perturba- 
tion strength increases. From now on, we define x as the 
partial sum <|I7} with D c = 800. 

Let us now look at the parameter dependence of the quantity 
yW (A) for finite temperature, plotted in Fig. [4] Recall from 
Ref. IU1I1 that x W characterizes the variation of the eigen- 
states, and as shown in Eq. (0 is a function of all level spac- 
ings or gaps. The plot is jagged due to the strong contributions 
of individual level crossings or near misses. However, one 
may identify two inflection points. For small couplings there 
is a local minimum of X > while for progressively larger cou- 
plings a global maximum appears. Past the global maximum, 
X^ 1 ' decreases as level avoidance begins to appear. From then 
on, as the couplings grow further x^ 1 -* continues to decrease. 

Our explanation is the following. For vanishing coupling 
the nearly degenerate subsectors of the spectrum give a large 
contribution to x^ due to small arguments to Gt(E n — E m ) 
and significant couplings | (n\ V\m)\ between neighboring lev- 
els. However, a slight increase of the coupling lifts those de- 
generacies, taking the system to the regime where perturba- 
tion theory applies since level crossings are rare (see Fig. |2^), 
and consequently x^ 1 ' is reduced. Continuing to increase the 
coupling, the energy levels between subsectors now begin to 
cross as one enters the so-called n-mixing regime. With suf- 
ficiently large coupling the system finally enters the quantum 
chaotic regime, where couplings between neighboring levels 
now lead to avoided level crossings, infrequent level cross- 
ings, and a consequent decrease in x^ ■ Note that the pertur- 
bative regime, for example, persists for a wider range of cou- 
plings for low energies than for higher energies, so the system 
is never entirely in one regime or another. This can be seen by 
comparing the energies and couplings with the classical chaos 



parameter r] — Ae 2 . For the hydrogen atom, higher energies 
result in a greater susceptibility of the electron to the external 
magnetic field over the Coulomb potential. 

In addition to taking a finite-temperature Gibbs state p = 
e~@ H /Z, another natural choice for the density matrix p is 
to consider the limit T — > oo. In this case our density ma- 
trix becomes proportional to a projector P over the subspace 
spanned by those eigenstates having well-converged eigenval- 
ues. Acting on this well-converged subspace, the state will 
have the form p^ = I/D c . The result for Xoo is plotted 
as a function of A in Figs. [4] and [5] With this form of the 
state we see that the contrast between the various regimes is 
enhanced. Moreover, this state uniformly weights the contri- 
butions of many more levels than the low-temperature case, 
suggesting that variation of Xoo may more closely mirror the 
level spacing statistics. It turns out that increasing the tem- 
perature strictly increases the magnitude of X , so that the 
"infinite"-temperature case of p = 1/D C has the largest value. 
Moreover, the location of the initial local minimum for small 
couplings shifts toward smaller values as temperature is in- 
creased. This is due to the sampling of higher energy levels, 
for which smaller couplings are required to traverse the vari- 
ous regimes. 

We can thus conclude that, as has already been found for the 
Dicke model in [9], the x^ P art of the OFS correctly incor- 
porates information about the level spacing statistics for the 
case of the nonlinearly coupled oscillators. As shown in Fig. 
[5] modifying the time t does not qualitatively change the plot 

of Xp ■ Noting the definition of Gt (x) following equation 
(|5j, where lim.t_ +00 Gt(x) = S(x), a larger time enhances the 
contributions to the sum due to level crossings or quasi- level 
crossings. The growth of x*- 1 -* for larger times thus reflects 
its sensitivity to the presence (absence) of level crossings, and 
hence the transition from regularity to chaos. 



C. The x (2) term 

So far we have treated the part of the OFS, x^ (A), which 
reflects the variation of the eigenstates and depends on the 
level spacings. The other term, X^W' i s proportional to 
the variance of the diagonal elements of the perturbation 
V = d\H with respect to p in the energy basis. In Fig. |6]one 
sees that x rapidly decreases with coupling strength. The 
values of x^ 1 ' an d X m both me regular and chaotic regimes 
differ by about two orders of magnitude, and this is not qual- 
itatively affected by temperature change. This consideration 
allows us to conclude that for this model the dominant con- 
tribution to the overall OFS is given by the fluctuations of the 
diagonal part of the perturbation V . However, though the con- 
tribution from x^ in this cas e is relatively small, it possesses 
a richer functional dependence on the coupling. In particular, 
while x/ 2 ' decreases mono tonic ally in the studied parameter 
interval, x^ 1 ' reflects the transition through the three described 
regimes. 

In the special case where the diagonal part of the pertur- 
bation in the energy basis, Vd, is zero, this second term will 
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FIG. 4: (Color online) The first term of the OFS, x (1) (A), with K = 
120, D c - 800, and t = 100 for (a) small temperatures and (b) 
higher temperatures. 
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FIG. 5: (Color online) The first term of the OFS, x (1) (A), with K = 
120, D c = 800, and various high temperatures at (a) t = 200 (b) 
t = 400. 
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FIG. 6: (Color online) The second term of the OFS, x (A), taking 
K = 120, D c = 800, and t = 100 



clearly vanish [8] . Thus, for an experiment or theoretical treat- 
ment in which it is desired that xh be isolated, choosing a 

(2) 

perturbation of this form will eliminate \p ■ 

Our results show that the behavior of the total OFS as a 
function of the parameter that drives the transition from the 
regular regime to a chaotic one encodes the already-mentioned 
resilience of the quantum evolutions to small non-random per- 
turbations even for specific systems. Indeed, for our model the 
statistical distinguishability between neighboring evolutions 
decreases with A, showing that the resilience of the system to 
perturbations dramatically increases when it becomes chaotic. 



V. CONCLUSION 

In this work we have used an operator-geometric quantity, 
the operator fidelity susceptibility (OFS), to study the transi- 
tion from regularity to quantum chaos for a pair of coupled 
two-dimensional harmonic oscillators. This model is equiva- 
lent to a hydrogen atom in a uniform external magnetic field, 
one of the prototypical systems for which both the regimes of 
classical and quantum chaos have been well-documented the- 
oretically and experimentally. We have seen that, by comput- 
ing the state-dependent part of the OFS, denoted x*- 1 ' (A), as a 
function of the coupling strength A one may distinguish three 
regimes in parameter space: perturbative, n-mixing, and quan- 
tum chaotic. A local minimum for small coupling strength 
corresponds to the boundary between the perturbative and n- 
mixing regimes, while the global maximum coincides with 
the occurrence of level crossings typical of a regular regime. 
As the level spacing statistics transform from Poisson (maxi- 
mum likelihood of level crossing) to Wigner-Dyson (vanish- 
ing probability of level crossings), X W correspondingly 
decreases. We may therefore conclude that x^-'(A) both in- 
corporates the information relative to the statistics of the spac- 
ings between neighboring energy levels and the information- 
theoretic notion of distinguishability between quantum evolu- 
tions. Our analysis shows that both the distinct elements of the 
OFS and the OFS as a whole are therefore effective tools for 
discriminating the different characters, i.e. regular vs. chaotic, 
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of the various regimes of a system's parameter space. 

We would like to remark that while the OFS approach may 
not be as efficient computationally as other more direct tech- 
niques, i.e. the level spacing statistics, the main goal in this 
work is to see whether the notion of resilience to perturbation, 
quantified by the OFS, is a useful one in the context of this 
well-known model. In the future, it would be interesting to 
experiment with this approach on models where the connec- 
tion between classical chaos and energy level statistics breaks 
down, such as in the odd-parity sector of the lithium atom 
12111 . Indeed, does the ability of \p to detect the transition 
from regularity to chaoticity solely depend upon a correspond- 



ing change in the level spacing statistics, or can it serve as a 
measure for chaoticity in systems without this property? 
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